A NON-STANDARD EVOLUTION PROBLEM ARISING IN 
POPULATION GENETICS 



FABIO A. C. C. CHALUB AND MAX O. SOUZA 

Abstract. We study the evolution of the probabihty density of an asexual, one locus pop- 
ulation under natural selection and random evolution. This evolution is governed by a 
Fokker-Planck equation with degenerate coefficients on the boundaries, supplemented by a 
pair of conservation laws. It is readily shown that no classical or standard weak solution 
definition yields solvability of the problem. We provide an appropriate definition of weak 
solution for the problem, for which we show existence and uniqueness. The solution displays 
a very distinctive structure and, for large time, we show convergence to a unique stationary 
solution that turns out to be a singular measure supported at the endpoints. An exponential 
rate of convergence to this steady state is also proved. 



1. Introduction 

A classical problem in population genetics is to study the evolution of a mutant gene. A 
standard approach to this problem is to consider a finite size population and to define a 
discrete dynamics for the evolution of the probability density of such a population. Usually, 
such models are Markov chains, in which the only absorbing states are the two pure ones. 
Therefore, one expects, for large time, convergence to one of these two states and, depending 
on which state is achieved, one says that the mutation has been either fixed or lost. For large 
populations, it is natural to ask for a continuous model that approximates this evolution. In 
a number of different ways, one arrives at a Fokker-Plank equation that describes either the 
evolution of the probability density (the so-called forward Kolmogorov equation), or what 
is sometimes called the transient fixation probability (the backward Kolmogorov equation). 
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From a mathematical point of view, it is interesting to notice that, for the fixation proba- 
bilities, it is easy to specify the appropriate non-homogeneous boundary conditions which, 
after subtraction of an appropriate multiple of a stationary solution, are recast as Dirichlet 
conditions. Nevertheless, this does not seem to be the case for the probability evolution. 
Since it must conserve mass, in many cases a condition of null probability current at the end- 
points is used (e.g. [22]). For a thorough introduction to the several aspects of mathematical 
population genetics, we refer the reader to the monographs by [3", 12] 

For a class of problems, however, these Fokker-Plank equations turn out to have degenerate 
coefficients at the boundaries, the classical Kimura equation (cf. [TYj) being the archetypal 
example. For the backward one, this is not a problem since the infinitesimal generator is, very 
generally, self-adjoint. For the forward equation, however, the underlying spectral problem is 
of the limit-point type and, thus, no boundary conditions can be enforced. In particular, one 
cannot control the flux of the solutions across the boundary of the corresponding domain, and 
the existence of conservation laws are not to be expected in general. This is an old issue in 
the study of diffusions, and it has been tackled by |13j . where the so called lateral conditions 
are derived, in order to ensure that the forward and backward equations are adjoint to each 
other. With these lateral conditions, however, the forward equation looses its differential 
character, and this led to a prevalence of the backward equation in the study of diffusions 
(particularly after [14j). We shall see below that is possible to ensure the duality of the 
backward and forward equations, while maintaining the differential character of the forward 
equation, within the framework of weak solutions. 

We shall study the forward Kolmogorov equation 



with F positive in (0,1), but with simple zeros at the endpoints, and with G vanishing 



F{x) = x{l — x), G{x) = x{l — x){r]x + /3) (forward Kimura with frequency selection; see [5]). 
Equation ([1]) is supplemented by the following conservation laws: 




dtpit, x) = dl (F(x)p(t, X)) - {G{x)p{t, x)) 
p(0, x) = p^{x) 



xG(0, 1), t>0 




(forward Kimura) and 



(2b) 



(2a) 




where ip satisfies 



(3) 



F{x)^" + G{x)^' = 0, 



V'(O) = 0, and V(l) = 1- 



More precise statements on the hypothesis made upon F and G are deferred to section [21 
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Remark 1. In population genetics, the function ^ is referred to as the fixation probability. 
Condition ()2ap is usually stated (or assumed), in the literature of population genetics, but 
condition ()2bp is not. These conditions have been derived in [5\, when obtaining the forward 
Kimura equation, with frequency selection, as a large population limit of the so called Moran 
process (cf. [23j/ 5*66 also [28j for an alternative approach. 

Before we proceed, we want to clarify the nature of the conservation laws given by ([2]). 
The backward equation and (formal) adjoint of ([1]) is given by 

{ dtf = F{x)dlf{t,x) + G{x)d^f{t,x) xG(0,l), t>0 
1 /(0,x) = /0(x). 

It is readily seen that any stationary solution to (jl]) is a linear combination of a constant 
and %l){x). Therefore, the conservation laws ([2]) are related to the kernel of the infinitesimal 
generator of dH). Finally, it should be mentioned that, if ([1]) is a correct approximation of the 
biological process, then one expects that the probability mass accumulates at the endpoints, 
as t goes to infinity [17] . 

The goal of this work is to clarify in what sense a solution to ([I]) that satisfies ([2]) exists, 
and how it behaves for large time. In contradistinction with [13] , which uses classical function 
spaces and has to modify equations ([T]) and ^ in order to obtain the duality relation, we shall 
always work with these equations, but in more general, non-normed, distributional spaces. 
This also differs from recent work in degenerate equations, as for instance: the controlability of 
degenerate heat equations [21j, with solutions in weighted Sobolev spaces; entropy solutions 
of Fokker-Planck from multilane traffic flow [10], where the conditions that might lead to 
concentration at the end points are explicitly avoided; and from the qualitative studies by 
[1]. See [15] for a general discussion of degenerate diffusion equations. We mention also the 
classical monographs [U |9] . 

Equation ([1]), with F{x) = x(l — x) and G{x) = x{l — x){r]x + (3) has been studied in 
reference [5], where a proof of existence and uniqueness in the sense of definition [1] is given, 
under the assumption of interior regularity. An announcement that also includes other results 
was made in [B]. See also [7j. More recently, the same problem has been studied through 
skillful, but formal, calculations in [22], with conditions of null probability current (formally) 
imposed. Thus, this work can be seen as complementary to the work by [13j by giving a 
differential formulation to the forward-backward duality for degenerate diffusions. Also, it 
can be regarded as an extension of [SJ [6], and as a rigorous proof of the formal calculations 
in [22]. 

The main results of the paper can be outlined as follows: let BM.^{[0, 1]) denote the space 
of (positive) Radon measures on [0, 1] . Then we have 
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Theorem 1 (outline). For a given E SA1^([0, 1]), there exists a unique solution p to 
Equation (QP, in a sense to be made precise in definition\^ with p e L~ ([0, oo); SA^+([0, 1])) 
and such that p satisfies the conservations laws The solution can he written as 

p{t, x) = q{t, x) + a{t)6o + b{t)6i, 

where 5y denotes the singular measure supported at y, and q G C°° (M^; C°°([0, 1])) is a 
classical solution to ([2]j. We also have that a{t) and b{t), belong to C([0, oo)) n C°°(M''"). In 
particular, we have that 

p e C~(M+;SX+([0,1])) n(:7°°(M+;C°°((0,l)) . 

For large time, we have that linit^oo (l{t^ x) = 0, uniformly, and that a{t) and b{t) are mono- 
tonically increasing functions such that: 

a°° := lim a(t) = / {\ — 'ip{x))p^ {x) dx and 

b°° := lim b{t) = [ 'il^{x)p^{x)dx, 

Moreover, we have that 

lim pit,-) = a°°<5o + 6°°5i, 
with respect to the Radon metric. Finally, the convergence rate is exponential. 

Remark 2. The coefficients of the singular measures, a{t) and b{t) are, respectively, the 
extinction and the fixation probabilities. Also, notice that the decomposition of p does not 
follows immediately from the linearity of ([T]j. As a matter of fact, neither of the summands 
are, per se, a solution to ([7]j in the sense of definition [7]. Heuristically, as ([ip is uniformly 
parabolic in each proper compact set of the unit interval, the parabolic operator erodes the 
interior density of the initial measure, which is then transferred into the boundaries and 
absorbed by the singular measures there. 

The outline of the paper is as follows: in section [2l we present background results for the 
classical (in a broad sense) solutions to ([T|). In section^ we introduce an appropriate definition 
of a weak solution and show that any solution of this type must satisfy the conservation laws 
We also show that, with this formulation, ([5]) is indeed the adjoint of ([I]). In section [U 
we present the proofs of existence and uniqueness. Section [5] discusses the convergence to the 
measures supported at the endpoints as time goes to infinity. 
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2. Preliminaries 

For the convenience of the reader, we present in this section some material that wiU be 
useful in the sequel. 

Let F, G : [0, 1] — > M be smooth, and assume that 

(1) F has single zeros at x = and at x = 1, and F{x) > 0, for x S (0, 1); 

(2) G has zeros at x = and x = 1. 

Hadamard's lemma (cf. [2j) then yields 

F{x) = x{l - x)^{x), ^'(x) > for X G [0, 1] and G(x) = x(l - x)n(x) 
Let us write, 

^(x) 



^'(x 

Then we can rewrite (HD as 



dtf = x{l-x)^{x)[d^f(t,x)+E{x)dJ{t,x)] XG(0,1), t>0 
/(0,x)=/0(x) 



(5) 

The stationary solutions of ([5]) are linear combinations of a constant and 

4,{x) = c-^ r Q-K^i')'^'- ds, c= [\-fo^^^)'i^ds. 
Jo Jo 

Existence of classical solutions to ([T|) can be established by Fourier series, and this is easier 

done by writing ([1]) in selfadjoint form. Let 

(6) e5r^(^)d^'u; = x(l -x)^'(x)p. 
Then ([T]) becomes 

(7) dtw = x(l - x)^'(x) ja^w - ^ [2E' + u;| . 

Remark 3. Since the standard maximum principle holds for C^'^ solutions of ([^, we find 
that, if the initial condition is nonnegative, then w{t, ■) is also nonnegative. Moreover, this 
holds also for p{t, •). 

Consider the associated spectral problem: 

+ I [2E' + E^]^ = Xeix)^, 

(8) 

Sturm-Liouville theory for singular problems allows us to conclude that ([8]) is a self-adjoint 
operator in ([0, l],0(x)dx), with a complete set of eigenfunctions. In what follows, all 
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spaces will be with respect to 9{x)dx, and we shall write (•, •) and || • II2 for the corrresponding 
inner product and norm, respectively. We also recall, see [U [25] for instance, that 

A. 

lim -4 = K. 

An important property of ([8]), which is proved in Appendix lA. I} is given by 

Lemma 1. The operator defined by |^ is positive-definite. 

We shall write (pj, j = 0, 1, 2, . . . , for the eigenfunctions of ([8]), with corresponding eigen- 
value Xj, and normalization ||(/5j||2 = 1- Also, for the spectral problem that arises in the 
original problem, we shall write 

(9) e^o-^'^'^'ipj = x{l - 

We shall need the asymptotic behavior of the eigenfunctions for large Xj. 
Lemma 2. There exists positive constants Ci and C2, independent of j, such that 

(10) W^jWoc < Ci and \\qj\\oo < C2Xf'^. 

The proof of Lemma [2] is given in Appendix IA.2I 

Finally, as in [27] for instance, we shall denote, for s > 0, the spaces 

00 I 

G L2([0,l],Mx)|^<KJ)AfVi e L2{[0,l],edx) \, ^{J) = 
j=o J 




j=0 

Since Radon measures are distributions of order less or equal to zero, we have — cf. [23 with 
minor modifications — that: 

Proposition 1. The initial value problem defined by Equation ^ and vu[0,x) = w^{x), with 
G ,BA4+((0, 1)) has the solution 

(11) w{t,x) = Y,^0{j)e-'^^ipj{x), w^ij) = {w'',ipj), 

which is unique in the class C°° (IR+; C°°([0, 1])) . 



Remark 4. It can be shown that, any standard weak solution definition to ^ will lead to 
the solution above — see for instance [HI [18]. Therefore, none of the conservation laws ^ 
can hold, and no classical-weak solution to |2H^ exists. 
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3. Weak solution and duality formulation 
We now make precise what we mean by a weak solution to ([T]). 

Definition 1. A weak solution to ^^ will be a function in L°° ([0, oo); SA^([0, 1])) that sat- 
isfies 

f oo />1 

p{t, x)dt(l){t, x)dxdt 



Jo 



CO /•! 

p{t, x)x{l - x)*(x) [dl(j){t, x) + E{x)d^(p{t, x)] dx dt 



^0 

1 





+ / p^{x)(f>{0,x)dx, 







where 



(j){t,x) eT = C^{[0,oo) X [0,1]) . 



Remark 5. Notice that the test functions in definition{l\ are required to be of compact support 
in [0, 1] and not in (0, 1) as usual. Similar definitions have been given in other contexts; see 
|19i I20j . where they are termed houndary-coupled weak solutions. 

Definition [I] can be recast in the framework of usual distribution theory, by identifying a 
Radon measure with a compactly supported distribution of nonpositive order (see [26lj- In 
this case, the distribution can act in C°°(M), but it is entirely determined by its behavior in 
the support; see for instance |16j . 

A glance at Definition [1] shows that, on the integral on the right hand side, the test function 
4> is applied to the operator on the right hand side of Thus, one could expect that any 
solution that satisfies ([T|) in the sense defined above, also satisfies the conservation laws ([2]). 

Proposition 2. Let p G L°° ([0, oo); BM-{[0, 1])). Ifx{x) is a stationary solution of then 
the quantity 

ri{t) = / x{x)p{t,x)dx 
Jo 

is constant in time. 

Proof. Let C(t) G C^((0,oo)). Then, (l){t,x) = C{i)x{x) is an appropriate test function. On 
substituting (j){t, x) in Definition [U we find that 

oo 

vit)c'it)dt = 0. 

Thus 7]{t) is constant almost everywhere. □ 
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Remark 6. We observe that standard spectral theory shows that both the infinitesimal gen- 
erators of ([1]) and dH) can be appropriately defined in a domain dense in L'^{{0, 1)) such that 
they are adjoints of each other. However, in this case, equation ([T]) will not be the forward 
Kolmogorov equation associated to On the other hand, in the sense of the pairing used in 
definitionUl gP with f{t, •) G C^{[0, 1]) is the adjoint of ^^ with p{t, •) G BM{[0, 1]). Thus, 
we recover the usual interpretation of the conservation laws given by the kernel of the adjoint. 



In what follows, it will be convenient to decompose a compact distribution, or a Radon 
measure, as the sum of a distribution without singular support at the endpoints, and two 
distributions singularly supported at the endpoints. We shall write £' to denote the space of 
compactly supported distributions in R. 

Lemma 3. Let v G E' , with supp(z^) = [0,1]. Then, the setwise decomposition 



where Vi is a compact distribution supported at x = i, and sing supp{p) C (0, 1). Moreover, 
if V is a Radon measure, then /u G B^A{{Q,1)), and Vi = Ci6i , with Cj G M, are singular 
measures with support at x = i. 



Proof. Let Q, i = 0,1,2 be a partition of unity in [0,1], subordinated to the open cover 
{[0, 2e), (1 - 2e, 1], (e, 1 - e)}. Let (p G C~([0, 1]). Define i^i, i = 0, 1 and by 



Then clearly v = VQ+fi+i^i. Also, it is readily seen that sing supp(//) C (0, 1). Since Co(^) = 1 
and Q^^'^x) = 0, n > 1, for x G [0,e), we find that uq is supported at x = 0, with a similar 
argument holding for z^i. Moreover, since a Radon measure is inner regular, the restriction 
of to (0,1) yields a Radon measure in (0,1). Finally, a Radon measure supported in a 



4. Existence and uniqueness 



[0,1] = {0}U(0,1)U{1} 



yields a decomposition of v, namely 



u = Uo + fl + Ui 




= 0,1 



singleton must be an atomic measure. 



□ 



For the initial condition, we shall write 




+ b^6 
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to denote the corresponding decomposition. Also, in order to show the existence of a solution 
to ([T|) in the sense of definition [H we shall temporarily consider p € -L°° ([0, oo); £"'), with 
support in [0, 1]. We shall write 

P = Po + q + Pi, 

for the decomposition of p. 

We now show that q must be, as a matter of fact, much more regular. 

Proposition 3 (Interior regularity). Assume that q° G eA1+((0,l)). // a solution to {iP 
exists, then q{t, x) must be the unique classical solution in the sense of section\^ with (7(0, x) = 
q^{x), i.e., 

00 

(12) 9(i,x) = 5;fo(i)(Z,e-^^*, 

3=0 

with qj given by Q, and q^{j) is j-th Fourier coefficient of q^. In particular, 

(7GC°°(E+;C°°([0,1])) . 
Proof Let </> G ([0, 00) x (0, 1)). Applying to definition [H we find 

fOO /"I 

q(t, x)dt(j){t, x) dx dt 



^0 



00 nl 



q{t, x)x{l - x)'^{x) [dl(j){t, x) + E{x)dx(l)it, x)] dx dt 
+ I g°(x)0(O,x)dx. 



JO 

00 





The result now follows by taking testing functions of the form 

^{t,x)=C{t)e~yo^^-)'i-4>{x), (eCc{[0,(x^)) and <^ G Ce((0, 1)), 

and then we use a standard Galerkin approximation procedure. □ 

Before we proceed, we observe that, since po and pi are distributions supported on a 
singleton, we must have, for some integers M and M' , that 

M M' 

(13) p{t, x) = Y, afe(t)4') + bk{t)5f^ + q{t, x), 

k=0 k=0 

(k) 

where (5^0 denotes the k-th distributional derivative of the singleton supported measure. 

Theorem 2 (Existence and uniqueness). The unique solution of HP in the sense of defini- 
tionUl with initial condition p^ G BA4^{[0, 1]) is given by 

p{t, x) = q{t, x) + a{t)5o + b{t)6i, 
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with q{t,x) given by (W(l . Moreover, we have 

a{t) = ^{{)) I g(s,0)ds + a° and h{t) = ^{l)f q{s,l)ds + }P . 
Jo Jo 

Proof. First, we define 

f={</>€C7~((0,«)) X [0,1])}. 

For Z > 0, we also define, 

%o = {0 G Cci[0, oo) X [0, md^Ht, 0) = 0, < n < /} , 

with a similar definition for Ti i. Notice that, for r > s, %fi C 

On substituting ([13]) in definition [H with cj) G T, using that q is smooth for t > and 
integrating by parts we obtain that 



M 



M' 



[(7(t,0)</.(t,0) +(z(t,l)</.(t,l)] dt+ 



^ ak{t)dtd',(t>{t, 0)+Y, hk{t)dtd^(t>{t, 1) 
+ / 5^afc(t)a^ {x{l-x)^{x)dl<l){t,x))\^^^ dt+ / ^afc(t)9^ (x(l - x)n(x)9,(/)(t, x))Uo dt+ 

+ / y^hk{t)dl {X{l-X)^{x)dlct>{t,x))l^^ dt+ / y^hu{t)dl (X(l-X)n(x)a,0(t,x))|,^, dt. 
JO Jo 



M 



M 



Restricting somewhat further, for e TM+ifi-, we find that 







aM{t)d,[x{l - x)^'(x)],=o9f +V(i,0) dt. 



Thus aM{t) = 0, and the sum can be only up to M — 1. Repeating the argument inductively 
yields M = 0. An analogous argument yields M' = 0. Thus only ao{t) and bo{t) can be 
nonzero. We now drop the subscripts and determine their values. Applying definition [1] to 
(p £ T, such that cpi^t, 1) = 0, we find that 

fco poo 

- a{t)dt(j){t,0)dt = g(t,0)^'(0)(/>(t,0)dt. 
Jo Jo 

Integrating by parts the corresponding relation for a{t), we obtain 



a{t)dt(l){t,0)dt 



l'(O) ^ qis, 0) ds + a° j dtcpit, 0) dt. 



Hence 



a{t) - ^{0) [ g(s,0)d, 
Jo 



s — = const, 
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everywhere, in as much as the integral is continuous. Since a(0) = a^, the identity follows. A 
similar calculation also shows that 



b{t) = ^{l) [\{s,l)ds + b°. 
Jo 



Uniqueness follows from proposition [3] and from the expressions for a{t) and b{t). Finally, 
notice that, since q{t,x) > 0, we have that both a and b are increasing. □ 

5. Large time behavior 

We now present some results for the behavior of the solution in the large time limit. 
Let us define 

6~ := [ il>{x)p^{x)dx = b^ + [ ip{x)q^{x)dx, 
Jo Jo 

-1 r-l 

p\x) dx - b°° = a° + / {l-iP(x))q°{x)dx. 
10 Jo 

Using the conservation laws ([2]), we also have 

6°°:= / il){x)p{t,x) dx = b{t) + / ■ilj{x)q{t,x) dx, 
Jo Jo 

p{t,x)dx -b'^ = a{t) + [ (1 - V'(x))g(t,x)dx. 
/o Jo 

Since < ip < 1 and q{t, •) > 0, we have that both a°° — a{t) and b°° — b{t) are nonneg- 

ative. From the representation given by p!T|) . we have that limt^oo |k(i,a^)||oo = 0. Hence, 

limt^oo CL{t) = a°° and lim^^oo b{t) = 6°°. 

Moreover, since q{t,x) > 0, we have 

a~ - o(t) + 6°° - = / q{t,x)dx=\\q{t,-)\\i, 

Jo 

which also yields the inequalities 

a°° -a(t) < ||g(t,-)lli and b°^ - b{t) < \\q{t, 
The behavior of the norm of q is given by the following result: 

Proposition 4. Let p be the solution to ^^ with an initial condition with q^ G BAi^{{0, 1)) 
and let Xq be the smallest eigenvalue of (E^. Then we have that 

lim e^«*||g(t,-)||i = Coo. 

In addition, if we assume that 

= x(l -x)*(x)e-3r=Wdy ^ ^7W+((0,1)) nP„ s > 0, 
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then there exists Cq^s > such that 

\\q{t,-)\\i<Co,s\\w^\\se-^^''. 
In particular, the same limit property and bounds apply to a°° — a{t) and b°° — b(t). 

Proof. For the first part, recall that 

oo 
j=0 

Let us write 

Qj = / qj{x) dx. 
Jo 

We observe that 

Q^. = (e-|/o^=(-)d^(^,)• 
However, since ez we do not have an immediate bound for \Qj\. 

On the other hand, we observe that qj satisfies 

-Xjqj{x) = dl [x{l - x)'^{x)qj{x)] - [x{l - x)Il{x)qj{x)] , 

which integrated yields 

*(0)g,(0) + M/(l)<7,(l) 



Q 



A, 



For large j, (jlOp guarantees that we must then have 



|Q,| <C2A7'/^ 



Thus, for t > 0, we have 



e"°1g(t,-)l|i = Qo^^°(0) + J]Q,^ro(j)e-(^^-^o)*, 

i=i 

and the result follows with Coo = Qqw'^{^)- 

For the second part, let us define the auxiliary functions: 

oo oo 

as{x) =^w^{j)\''J'^ipj{x) and ()s{t,x) = ^QjX^''''^ Lpj{x)e 

3=0 3=0 

Then, we have that 

\\q{t,-)\\i = {as,(3s{t,-)) < ||a,||2||/?.(t,-)ll2 = 

= ||t/;°||.e-^°*||e^o*/?,(t,-)||2 < Co,.||u;0|Ue-^°*, 

withC7o,« = Il/3s(0,-)ll2. 



□ 
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Theorem 3 (Exponential convergence). Let p denote the Radon metric, and let 

= a°°(5o + 6~5i. 
Under the same hypothesis of proposition^ we have that 

(14) lim e^«V(P,P°") < 2Coo. 

t— >oo 

With the additional hypothesis, we have that 

(15) p{p,p'^)<2Co,s\\uj^\\se-^''. 

In particular, Ji^D implies convergence in the Wasserstein metric. 

Proof. Recall that 



p{u,p)=sup< / f{x)d{u-p) 







/GC([0,1];[-1,1]) . 



But, for such / we have that, when t > 0. 



/(x)d(p~-p(t,x)) 







< / \d{p^-p{t,x))\ 
Jo 

< [ (a°° -a(t))5odx+ [ (6°° - 6(t))5idx + [ \qit,x)\dx 
Jo Jo Jo 

= a°° -a{t)+b°^ -b{t) + \\q{t,-)\\i 
= 2|k(t,-)lli- 

Now, both (jl4p and (jlSp follows from Proposition 4. □ 



Remark 7. In many applications, the slowest decaying mode ipQ is taken to he a quasi- 
stationary distribution for the diffusion process. The constant Coo is then the total probability 
mass of such a distribution. 

Appendix A. Postponed proofs 
A.l. Proof of positive-definiteness of ([8|). 

Proof of Lemma [II Let 

(16) i; = e-5r=(^)d^(^. 

Then ^ becomes 

(17) -v" -Ev' = Xe{x)v, t;(0) = w(l) = 0. 

When A = 0, then (jl7p becomes the stationary version of ([5]), with Dirichlet boundary 
condition. Its general solution, v, is given by D = ci + C2ip, which does not satisfy the 
required boundary conditions. Thus, zero cannot be an eigenvalue of p!7|) . Moreover, since 
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the transformation (jl6p preserves the oscillation properties of the eigenf unctions, we have 
that the eigenfunction vq, corresponding to the smallest eigenvalue Aq, will not have any 
zeros inside (0, 1). Let us assume, without loss of generality, that > in (Oi !)• It must 
have a point of maximum x = x* G (0, 1), where Vq{x*) = 0. Hence we must have 

vojx*) „ 
x*[l — x*)^f[x*) 

Note that Vq{x*) ^ 0, otherwise we would have Aq = 0. Since it is a maximum, we must have 
Vq(x*) < 0. Since, vo{x*) > 0, we have Aq > 0. □ 

A. 2. Proof of the asymptotic estimates. 

Proof of Lemma \M For the proof, we drawn on results by [231 chapter 12] that are summarized 
as follows 

Theorem 4. Let 

Also let 

'1/2 



where Ji is the standard Bessel function of order one, and Aqj is choosen such that ||<^j||2 = 1- 
For large j, we have that 

- ^,|U < i^iFi(Af [Cl^/^)exp [pF2{C?j 

where K\, K2 are positive constants andFi, F2 are positive and hounded continuous functions. 

With this result, we can now prove the asymptotic behavior for ipj and qj 
Let 

and Ji is the standard Bessel function of order one. Let z = Then we find that 

where = Aj''^|C(l)j- For large j, we have from the asymptotic behavior of Ji at infinity 
that 

= C\]'^ + 0(1). 
Also, since uJi{u) < Cy^u), for large u, we have 
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Combining these two results, with Theorem SI we have the first result in ([10 
For the second one, we observe that 



dcy/' Ji(£) 

dx J z 



Combining with the estimate for Aqj, we have the result. □ 
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